Energy landscape, two-level systems and entropy barriers in Lennard-Jones clusters 
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We develop an efficient numerical algorithm for the identification of a large number of saddle points 
of the potential energy function of Lennard- Jones clusters. Knowledge of the saddle points allows 
us to find many thousand adjacent minima of clusters containing up to 80 argon atoms and to locate 
many pairs of minima with the right characteristics to form two- level systems (TLS). The true TLS 
are singled out by calculating the ground-state tunneling splitting. The entropic contribution to all 
barriers is evaluated and discussed. 



The use of an effective potential energy to describe the 
dynamical properties of disordered systems has become a 
very powerful tool of investigation. In particluar, there is 
broad consensus on the idea that important relaxational 
processes can be described by transitions between ad- 
jacent potential-energy minima in the multidimensional 
configuration space. The latter, however, is so complex 
even for systems containing a moderate number of atoms, 
that singling out relevant pairs of minima has required 
important numerical efforts and skill by many authors 
. One of the goals of this kind of research is the iden- 
tification of pairs of minima suitable to produce two-level 
systems (TLS); these are the pairs of nearly degenerate 
minima separated by a barrier that, following tunneling, 
produces a total splitting of the order of «< 1 K, thus 
explaining the anomalous low-temperature specific heat 
and thermal conductivity of glasses ||-[l(J . 

Identifying adjacent minima is not simple by molecular 
dynamics methods; it is clear that the key stationary 
points arc not the minima but rather the (first-order) 
saddle points, descending from which adjacent minima 
are automatically found. But there are other reasons why 
the search for saddle points looks desirable, which are in 
a way or another related to the fact that the minimum- 
minimum relaxation process takes place along a classical 
path that passes close to the saddle, i.e. the least-action 
path. 

The first reason is that evaluating the TLS ground 
state splitting is a quantum mechanical problem, but the 
solution of the Schroedinger equation in 3iV dimensions 
is out of question and one is forced to think of ways of 
drastically reducing the dimension of the problem, possi- 
bly to one. It is therefore necessary to have an indication 
about how much and under which circumstances this re- 
duction is feasible; as we will discuss later this kind of 
information requires a detailed knowledge of the path 
along which tunneling is considered to take place and of 
the potential energy along it. 

The second reason, which is related to the first one, 
is that relaxation processes of all kinds, either the 
low-temperature ones that produce TLS or the high- 



temperature ones that govern the glass transition, in 
principle are controlled not only by energy barriers, but 
also by entropy barriers, i.e. by the topological difficulty 
of proceeding from one minimum to another through a 
narrow route even in the absence of potential energy bar- 
riers |Ii}-[l| . Note that entropy barriers would also affect 
the quantum mechanical splitting of the ground state in 
much the same way as energy barriers do fli|] . 

Heuer has recently studied the potential-energy 
topology of systems containing 32 atoms placed in a 
box with periodic boundary conditions and interacting 
through a Lennard- Jones- like potential introduced in ref. 
jjj. At a numerical density p = 1 in units of the nearest 
neighbor distance, the author finds 367 minima of dif- 
ferent energy, which become 75 by decreasing the box 
size by ~2%. These values are exceedingly small when 
compared to the « 10 11 different minima expected for 
a Lennard-Jones cluster of 32 atoms jl5|]: the bound- 
ary conditions have a dramatic effect on the number of 
minima (and probably of all kinds of stationary points). 
Heuer [Q notes that performing the simulation at vari- 
able box length produces many more minima; in view 
of this situation, we think it is very important to search 
for the stationary points in the absence of any boundary 
condition, i.e. for free clusters. The reason is that even 
periodic boundary conditions with variable density are 
likely to introduce spurious correlations among different 
boxes, that might affect the results. The price to be paid 
is the presence of surface effects which (like the correla- 
tions of the boundary conditions case) become less and 
less important as the system size is increased. 

The systems we studied are Ar clusters with up to 
N =80 atoms, interacting through the Lennard-Jones 
potential 

% j>i J J 

with a = 3.405 A and e/A: B =I25.2 K. 

Since we think of our clusters as parts of a glass, in 
which case translations and rotations are forbidden, we 
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decided to set six coordinates to zero (one atom is fixed 
in the origin, one on the x axis and one on the xy plane) 
in order to eliminate translations and rotations. Free 
translations of the system as a whole are decoupled by 
other degrees of freedom, but in non-rigid bodies like 
clusters rotations and vibrations are coupled. Therefore, 
the above conditions leave the potential energy of the 
stationary configurations unaffected because it depends 
only on mutual distances, but the elimination of rotations 
affects vibrational frequencies by altering the effective 
mass. We do not expect any systematic bias of the results 
because the coordinates to be fixed are chosen at random. 

The saddle points were found with the procedure de- 
scribed in detail in ref. , and which consists of the fol- 
lowing steps: (i) descent towards a minimum by the con- 
jugate gradient method, starting from a randomly chosen 
configuration; (u) ascent towards the vicinity of a saddle 
following the eigenvector corresponding to the minimum 
eigenvalue of the Hessian; this requires diagonalization of 
the Hessian at each step; (Hi) once the potential energy 
along the path of the previous item starts decreasing, we 
take the corresponding configuration as the starting point 
for a Newton- Raphson stationary point search |l6| . This 
means that the non linear system of equations 
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is solved by successive iterations for i=l...(3iV-6). The 
steps (i — ii) are required in order that the Newton pro- 
cedure converges; the iteration ends up in a first order 
saddle in about 30% of cases and the saddle is located 
with extreme accuracy. The next step is a minimum- 
eigenvalue descent on both sides of the saddle, which 
provides approximate adjacent minima that are subse- 
quently fed into the Newton algorithm for accurate lo- 
cation. With this numerical procedure we found many 
thousand minimum-saddle-minimum triplets for several 
values of N in the range 6-80. The number of triplets 
found and of those with asymmetry Ajl K, are reported 
in the second and third columns of Table I, respectively. 
The last 2 columns of the Table entail the evaluation of 
the ground state splitting, which will be discussed later. 

The first check we performed was to look for corre- 
lations between the asymmetry and the barrier height; 
contrary to previously reported results on periodic sys- 
tems [^,^1 we found no correlation. The pairs of minima 
were further characterized by their euclidean distance, d, 
and by the participation number, N p = d /d^uxi where 
dmax is the displacement of the atom that moves most. 
The distribution of d is peaked at w 4cr for N — 13 — 42 
and at « 6<r for N — 80, with tails extending up to 
« 15cr. These tails reflect the free- surface effects that 
we know affect the clusters, but interestingly enough the 
TLS pairs have in all cases d <« 2a. This is convenient 
for an a priori selection of candidate TLS. The partici- 
pation number for N = 42,80 is reported in Fig. 1. In 



column (b) of Table I we report the number of candidate 
TLS, i.e. those triplets that have both Ajl K and the 
barrier higher than the lowest vibrational eigenvalue of 
either minimum; the distribution of their participation 
number is not very different from that of all pairs for 
N < 42, while for N — 80 we have too few pairs to draw 
any definite conclusion. 

In order to identify the TLS, we need the total split- 
ting of the ground level, resulting from asymmetry, A, 
and tunneling, St- As mentioned above, the evaluation 
of the latter is a quantum mechanical issue and can be 
solved relatively easily only if the many- dimensional dy- 
namics can be reduced to a 1-dimensional one. If this 
is at least approximately feasible, then the computation- 
ally most convenient way of evaluating St is to use the 
semiclassical WKB approximation (ljjl^l that provides 
an explicit expression in terms of the potential energy 
profile V(x): 
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with D = [1 + exp(2S)] 1 and the action integral S = 
I / {ot[2F(x) — Tiuj^^dx is evaluated between the clas- 
sical turning points a and 6; u) is the classical frequency 
in one well. The total splitting is then obtained as 
S ~ y/ A 2 + S'ip. As mentioned, in principle it is possi- 
ble to use equation (2) for evaluating the tunneling split- 
ting only if the Schroedinger equation can be separated 
into independent ones for each degree of freedom p9| ; 
we shall assume that (2) can be employed under the less 
restrictive condition that the relevant classical degree of 
freedom, i.e. the coordinate along the least action path, is 



decoupled by the remaining ones. Demichelis et al J14 20 
first found that this actually occurs for argon with peri- 
odic boundary conditions, and the same has been veri- 
fied for candidate TLS in clusters pl| : the higher vibra- 
tional eigenvalues, obtained by diagonalizing the dynam- 
ical matrix at configurations along the least-action path, 
do not change appreciably along the path itself, while the 
ground one obviously does. Under these circumstances 
the use of Eq. (2) is reasonably justified [|19f| . The least- 
action path was determined by minimizing the action in- 
tegral S following the procedure of ref. 14 1, and by tak- 
ing as initial path the one starting from the saddle and 
reaching the minima by a minimum eigenvalue search. 
The action integral was evaluated also along the straight 
line connecting the minima; for the TLS the straight-line 
integral turns out to be only a factor 2-5 larger than the 
least-action one. 

In the last column of Table I (TLS) we report the num- 
ber of triplets that have a total (asymmetry plus tunnel- 
ing) splitting in the ground state S < 1 K, and thus have 
all the characteristics of TLS 22 1 . We see that interest- 
ingly enough the category of TLS almost coincides with 
that of column (6) . But most remarkably, it seems that 
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n(TLS)/n(Total) is roughly of the order of 10 -3 and ap- 
parently tends to increase with N; if this trend should 
be confirmed it would mean that surface effects do not 
favor TLS, so that the failure to observe them in periodic 
systems of comparable or larger p(J size should prob- 
ably be ascribed either to correlations or to density. In 
any case, visual inspection of the TLS relative to N = 80 
shows that surface effects are still very important because 
in all 3 cases there is at least one surface atom that moves 
most. For N = 42, 80 we verified that TLS are isolated, 
in the sense that we did not find any low-barrier route 
that connects one minimum of the pair to a third almost 
degenerate minimum. This was to be expected since only 
w 1/100-th of the pairs are nearly degenerate (Table I, 
column (a)). 

The entropic contribution to classical relaxation was 
considered along the following lines. The classical one- 
dimensional transition rate between minima separated by 
a barrier Eb is 
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where w is the vibration frequency in the potential well; 
in many dimensions and under some assumptions con- 
cerning thermodynamic equilibrium and the absence of 
back crossing, it is found that the entropic-barrier effect 
of degrees of freedom other than the single considered 
path can be accounted for by the following substitution 
in the pre-exponential factor of (3) pF 
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Here ojf 1 and uij are the vibrational eigenfrequencies at 
the minimum and at the saddle respectively, and the 
prime indicates that the imaginary frequency at the sad- 
dle point has to be omitted from the product. 

The entropic ratio R has been calculated for all 
minimum-saddle-minimum triplets for N — 42 and 80, 
and is reported for N = 80 in fig. 2 as a function of 
the barrier. For N =42 the distribution is qualitatively 
similar but much more scattered. As can be seen, the 
majority of the triplets have R in the range 1-10 -2 , with 
a high-barrier tail extending down to R rs 10~ 6 . Only 
relatively few triplets (0.45%) have R > 1 (but in any 
case not exceeding the value 10); for N —42 this fraction 
becomes 0.75%. The present results seem to indicate 
that entropic barriers are not especially active in slowing 
down the relaxation dynamics, and that the effect may 
become less important with increasing N. This result is 
contrary to what found recently for a different model of 
fragile glass with periodic boundary conditions pq] : the 
role of boundary conditions on the results of simulations 
on relatively small systems definitely deserves further in- 
vestigation. 

For the TLS, in agreement with the reported indepen- 
dence of the eigenvalues other than the ground one on 



the position along the least-action path [£0|] , we find R w 
1. Considering that for N —80 R is the ratio of two prod- 
ucts of 234 numbers each, the result R « 1 is remarkable 
and shows that TLS belong to a rather peculiar class of 
triplets. In this case, entropic barriers do not play a rel- 
evant role in determining the quantum splitting of the 
ground level, though in general they quench the splitting 
like ordinary potential barriers. 

In summary, TLS have been shown to exist in Lennard- 
Jones clusters and have been characterized; the order of 
magnitude of the ground level quantum splitting has been 
evaluated by means of the semiclassical 1-dimensional 
WKB formula; topological hindrances do not play a sig- 
nificant role in determining this splitting. Moreover, the 
relaxation dynamics of the vast majority of minimum- 
saddle-minimum triplets appears to be little affected by 
entropic barriers. As far as we are aware, this kind of 
information on tunneling and relaxation processes has 
not been provided before, and it was made possible by 
the use of a numerical procedure that produces a large 
amount of adjacent minima. The quantitative applica- 
tion of these results to glasses requires great caution, be- 
cause surface effects are certainly present in our clusters. 
However, the influence of boundary conditions and den- 
sity on the properties of the landscape |?J is also very 
great and the present zero-pressure and correlation-free 
results may certainly help in getting a closer insight into 
the problem. 

We are grateful to A. Ranfagni and D. Mugnai for very 
useful discussions. 
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TABLE I. Number of pairs of adjacent minima found for 
N = 6 — 80. A is the asymmetry and 5 the total splitting of 
the ground state (see text), (a): A < 1 K; (b): Same as (a) 
but with barrier higher than minimum vibrational eigenvalue; 
(TLS): same as (b) but with total splitting 5 < 1 K and 
tunneling splitting 8 T > 10" 15 K (see ref. [20]). 
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FIGURE CAPTIONS 



Fig. 1. Participation number for clusters with N = 42 
(squares) and 80 (circles). 



Fig. 2. Entropic ratio R, as defined in Eq. (4), for 
N = 80 as a function of the saddle-minimum energy bar- 
rier; full circles: upper minimum; open squares: lower 
minimum. 
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" Energy landscape, two level...." 

Figure 1 
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